详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

内装式空射运载火箭重力出舱机箭耦合动力学分析

研究了以运输机为平台的内装式运载火箭空射过程载机和火箭的耦合动力学建模。建模针对两个阶段：第一阶段，火箭固定于载机机舱内，两者构成一个整体，按照普通刚体的力学方法处理；第二阶段，火箭解锁后，沿着舱内的发射筒向外滑行，载机和火箭形成两刚体相互作用的耦合系统，基于牛顿-欧拉法建立系统动力学模型。载机在空射火箭过程中，油门和升降舵满偏，在加速前飞的同时拉大姿态俯仰角，火箭在自身重力分量和惯性力的作用下，沿着机舱内的发射筒加速向外滑行，直至与载机分离。数值仿真分析了空射过程载机的重要力学参数的变化过程，验证了载机操控策略的可行性和安全性，可为未来中国空射运载火箭技术设计提供数据参考。     

Effect of bristle pack position on the rotordynamic characteristics of brush-labyrinth seals at various operating conditions

Over the last few decades, the research on the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals is not sufficient. To this end, two kinds of brush-labyrinth seals for the bristle pack element installed upstream of the labyrinth teeth named BSU and installed downstream of the labyrinth teeth called BSD were used to investigate the effect of bristle pack position on the rotordynamic characteristics of the brush-labyrinth seals. Using the numerical model combining the porous medium model and the whirling rotor method, the rotordynamic characteristics of the BSU and BSD at various operating conditions including four kinds of pressure ratios, five kinds of inlet preswirl speeds and four kinds of rotor spinning speeds were conducted. The obtained results show that the effects of operating conditions on rotordynamic coefficients for the different seal configurations are different. The direct stiffness, cross-coupled stiffness and direct damping of the BSU are lower than those of the BSD. The rotordynamic coefficients of the BSU are more insensitive to the operating conditions variation. From the perspective of the seal stability, the BSU is a better brush-labyrinth seal configuration at high pressure ratio, high positive preswirl or high rotor spinning speed conditions. While in the case of low pressure ratio, low positive preswirl or low rotor spinning speed conditions, the BSD is a better choice.     

基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

A progressive approach to predict shot peening process parameters for forming integral panel of Al7050-T7451

In this paper, a progressive approach to predict the multiple shot peening process parameters for complex integral panel is proposed. Firstly, the invariable parameters in the forming process including shot size, mass flow, peening distance and peening angle are determined according to the empirical and machine type. Then, the optimal value of air pressure for the whole shot peening is selected by the experimental data. Finally, the feeding speed for every shot peening path is predicted by regression equation. The integral panel part with thickness from 2 mm to 5 mm and curvature radius from 3200 mm to 16000 mm is taken as a research object, and four experiments are conducted. In order to design specimens for acquiring the forming data, one experiment is conducted to compare the curvature radius of the plate and stringer-structural specimens, which were peened along the middle of the two stringers. The most striking finding of this experiment is that the outer shape error range is below 3.9%, so the plate specimens can be used in predicting feeding speed of the integral panel. The second experiment is performed and results show that when the coverage reaches the limit of 80%, the minimum feeding speed is 50 mm/s. By this feeding speed, the forming curvature radius of the specimens with different thickness from the third experiment is measured and compared with the research object, and the optimal air pressure is 0.15 MPa. Then, the plate specimens with thickness from 2 mm to 5 mm are peened in the fourth experiment, and the measured curvature radius data are used to calculate the feeding speed of different shot peening path by regressive analysis method. The algorithm is validated by forming a test part and the average deviation is 0.496 mm. It is shown that the approach can realize the forming of the integral panel precisely.     

机载托架自动校准系统的设计与实现

为了解决机载托架传统校准方法精度低、耗时、费力等问题,设计了一种高精度、高效率并且操作简单的成品托架自动校准系统。该系统运用调平精度较高的"循环多次"最高点不动调平方法,建立高精度机载托架校准的静力学数学模型,运用VC与Matlab的COM接口编程技术开发机载托架自动校准系统,并建立良好的人机交互界面。试验结果表明,利用该方法对机载托架进行校准,无论是校准精度还是校准效率都较传统方法有很大提高。     

Fault mode probability factor based fault-tolerant control for dissimilar redundant actuation system

This paper presents a Fault Mode Probability Factor (FMPF) based Fault-Tolerant Control (FTC) strategy for multiple faults of Dissimilar Redundant Actuation System (DRAS) composed of Hydraulic Actuator (HA) and Electro-Hydrostatic Actuator (EHA). The long-term service and severe working conditions can result in multiple gradual faults which can ultimately degrade the system performance, resulting in the system model drift into the fault state characterized with parameter uncertainty. The paper proposes to address this problem by using the historical statistics of the multiple gradual faults and the proposed FMPF to amend the system model with parameter uncertainty. To balance the system model precision and computation time, a Moving Window (MW) method is used to determine the applied historical statistics. The FMPF based FTC strategy is developed for the amended system model where the system estimation and Linear Quadratic Regulator (LQR) are updated at the end of system sampling period. The simulations of DRAS system subjected to multiple faults have been performed and the results indicate the effectiveness of the proposed approach.     

Topology optimization in lightweight design of a 3D-printed flapping-wing micro aerial vehicle

Topology optimization is an effective method to obtain a lightweight structure that meets the requirements of structural strength. Whether the optimization results meet the actual needs mainly depends on the accuracy of the material properties and the boundary conditions, especially for a tiny Flapping-wing Micro Aerial Vehicle (FMAV) transmission system manufactured by 3D printing. In this paper, experimental and numerical computation efforts were undertaken to gain a reliable topology optimization method for the bottom of the transmission system. First, the constitutive behavior of the ultraviolet (UV) curable resin used in fabrication was evaluated. Second, a numerical computation model describing further verified via experiments. Topology optimization modeling considering nonlinear factors, e.g. contact, friction and collision, was presented, and the optimization results were verified by both dynamic simulation and experiments. Finally, detailed discussions on different load cases and constraints were presented to clarify their effect on the optimization. Our methods and results presented in this paper may shed light on the lightweight design of a FMAV.     

单向电永磁作动器结构优化与动态特性分析

针对电磁作动器的使用要求，设计出一种具有能耗低、吸力大及响应快等特点的新型单向电永磁作动器。结合有限元仿真软件，利用有限元法研究作动器的结构参数如铁芯齿形结构、平面结构、无永磁结构、铁芯的倾斜角度及衔铁位移（气隙厚度）等因素对作动器电磁吸力和响应速度的影响。结果表明，电磁吸力随着作动衔铁位移的增加而大幅度减小；倾角结构的铁芯设计能使作动器性能达到最佳，且铁芯的倾斜角度为50°时，电磁吸力有最大值3071N。     

反推力装置运动学与动力学仿真

为了研究叶栅式反推力装置各部件在工作过程中的运动学与动力学特性,根据机构运动原理对反推力装置进行简化并建立了运动学与动力学数学模型。以滑动整流罩位移与阻流门所受气动负荷为输入进行运动学与动力学仿真,得到了反推力装置各部件的位移、速度及受力特性曲线,并对比分析了在不同尺寸参数下各部件特征点的运动轨迹和反推力装置负载力变化。结果表明:运动学及动力学仿真结果与工程实际相符；反推力装置机构参数选择不合理时,各设计点会发生干涉现象并导致机构无法运动;机构参数变化对负载力最大值影响尤为突出,在阻流门AC段长度值增大6%,阻流门CB段长度值减小9%的情况下,负载力正向最大值将增大19.53%,负向最大值增大12.67%。研究方法及研究结果可为反推力装置运动学及动力学分析,以及为反推力装置机构优化设计提供参考。